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We study the behavior of the dual quark condensate Ei in the Nambu-Jona-Lasinio (NJL) model 
and its nonlocal variant. In quantum chromodynmics Ei can be related to the breaking of the center 
symmetry and is therefore an (approximate) order parameter of confinement. The deconfinement 
transition is then signaled by a strong rise of Ei as a function of temperature. However, a similar 
behavior is also seen in the NJL model, which is known to have no confinement. Indeed, it was 
shown that in this model the rise of Ei is triggered by the chiral phase transition. In order to shed 
more light on this issue, we calculate Ei for several variants of the NJL model, some of which have 
been suggested to be confining. Switching between “confining” and “non-confining” models and 
parametrizations we find no qualitative difference in the behavior of Ei, namely, it always rises in 
the region of the chiral phase transition. We conclude that without having established a relation to 
the center symmetry in a given model, Ei should not blindly be regarded as an order parameter of 
confinement. 
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I. INTRODUCTION 


Spontaneous chiral symmetry breaking and confine¬ 
ment are among the most important features of nonper- 
turbative low-energy quantum chromodynamics (QCD). 
Both properties have been studied from first principles 
within lattice QCD, which nowadays provides an excel¬ 
lent description of the meson and baryon spectrum in 
vacuum [T]. Turning to nonzero temperature T, lattice 
QCD predicts a rapid but smooth crossover from the 
confined phase with broken chiral symmetry at low T 
to a deconfined phase with (approximately) restored chi¬ 
ral symmetry at high T m E]. Moreover, it is found 
that both transitions take place in the same temperature 
range. The corresponding order parameters are the chiral 
quark condensate (V'/) and the Polyakov-loop (PL) ex¬ 
pectation value HE], respectively, which strictly speak¬ 
ing belong to opposite limits of the theory. While (/’/) 
is a strict order parameter for chiral symmetry breaking 
in the limit of massless quarks, the PL expectation value 
is a strict order parameter for confinement in the limit of 
infinitely heavy quarks, where it can be related to Z(3) 
center-symmetry breaking and thereby to the free energy 
of a static quark-antiquark pair HII]. 

The identification of an appropriate order parameter 
for the confinement-deconfinement phase transition in 
QCD with finite (i.e., non-infinite) quark masses stands 
on less solid grounds. In this context the authors of 
Ref. [5] have proposed the dual quark condensate Ei as 
an alternative. Starting point is the generalized quark 
condensate ('i/'V')^, which is the analogue of the usual 
quark condensate, but evaluated for quark fields with 


twisted boundary conditions 

tp{x,(3) = e~"^tp{x,0) (1) 

in the imaginary time direction. Here /3 = 1/T is the 
inverse temperature and / G [0, 2tt) is an arbitrary an¬ 
gle. Hence, physical fermions, which obey antiperiodic 
boundary conditions, correspond to / = tt. 

The dual quark condensate E„ is then defined as the 
Fourier transform with respect to /, 

= ( 2 ) 

where n is an integer. In lattice representation this can 
be written as a sum of Wilson loops winding n times 
around the temporal boundary [Sj. In particular, since 
the PL is the shortest loop with winding number 1, the 
case n = 1 may be viewed as a collection of generalized 
Polyakov loops with spatial fluctuations. Therefore Ei 
has been termed “dressed Polyakov loop” [5]. 

El and the ordinary (“thin”) PL transform in the 
same manner under center transformations, which mo¬ 
tivates the consideration of Ei as an order parameter 
for the deconfinement phase transition. Moreover, since 
the spatial fluctuations are suppressed for infinite quark 
masses, Ei reduces to the thin PL in this limit. On the 
other hand, as seen from its definition, it is also related 
to the quark condensate, albeit with unphysical bound¬ 
ary conditions. This hints for a possible connection be¬ 
tween chiral and deconfinement phase transition, explain¬ 
ing why both transitions occur in the same temperature 
region 00- 

Another important feature of Ei is that it is not re¬ 
stricted to lattice formulations of QCD, but it can also be 
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calculated within continuum approaches, like the func¬ 
tional renormalization group method [9] or Schwinger- 
Dyson equations nnHH, where the calculation of the 
thin Polyakov loop is not possible in a straightforward 
manner. These investigations have confirmed that chi¬ 
ral restoration and deconfinement phase transitions take 
place in the same regime, even when the analysis is ex¬ 
tended to a nonvanishing chemical potential [Umi, a 
region which is not accessible in lattice QCD because of 
the sign problem. 

In addition to these studies directly rooted in QCD, 
Si has also been investigated within effective models of 
strong interactions. This was done first in the Polyakov- 
loop extended NJL (PNJL) model [1^1 and a bit later 
in the standard NJL model m- It was found in both 
cases that the qualitative behavior is similar to the QCD 
results, i.e., the rise of Si, which in QCD signals the 
onset of deconfinement, is found approximately in the 
same temperature region where chiral restoration takes 
place. This is particularly remarkable for the NJL model, 
which has been invented as a model for chiral symmetry 
breaking [niiiH] but is known to feature no confinement. 
Indeed, since there are no gauge fields and, hence, no 
center symmetry in this model, a connection between Si 
and confinement cannot be made{^ Instead, a Ginzburg- 
Landau type analysis revealed that in the NJL model 
the behavior of Si is triggered by the chiral phase tran¬ 
sition po] . 

The situation is somewhat less clear for the PNJL 
model, where confinement effects are included by cou¬ 
pling a gluon background to the quarks [mug. How¬ 
ever, in this way confinement is realized only statisti¬ 
cally, i.e., quark effects to thermodynamic quantities are 
suppressed at low temperature, but the quarks remain 
as physical states in the spectrum. As a consequence 
mesons can decay into quark-antiquark pairs |23j . which 
should not be possible in a confining theory. 

Already before the PNJL model was invented, there 
have been various attempts to model confinement by 
modifying the analytic properties of the quark propaga¬ 
tor, e.g. P3H1Z]- One possibility is a quark propagator 
without real singularities in the time-like momentum re¬ 
gion. Formally, this is related to the violation of reflec¬ 
tion positivity |281129) , meaning that quarks do not exist 
as physical states in the particle spectrum and are thus 
confined. 

As mentioned above, the NJL model in its original for¬ 
mulation does not support confinement; the structure of 
the quark propagator in this model is consistent with that 
of a free particle. However, employing the proper-time 


^ A similar conclusion was drawn in Ref. m, where Si has been 
explored in three-dimensional quantum electrodynamics, a con¬ 
fining theory which, however, does not possess a center symme¬ 
try. Therefore, although the behavior of Si is qualitatively the 
same as in QCD, it cannot directly be linked to the confinement 
transition. 


regularization, the unphysical quark production thresh¬ 
old can be avoided by introducing an infrared cutoff, as¬ 
sociated with the confining scale 1301 in]. The resulting 
propagator does neither develop real nor complex poles. 
This is another statement of confinement in the sense 
that the excitation described by a pole-less propagator 
can never reach its mass shell. 

Alternative attempts to simulate confinement utilize 
nonlocal extensions of the NJL model (nNJL model) [3^ 
l35] . Thereby the interaction is designed in such a way 
that the quark propagator has no real but complex con¬ 
jugate poles0 These complex singularities have been in¬ 
terpreted as confined quasiparticles [Ml EH- Through the 
incorporation of temperature into the model, the complex 
singularities may turn real and a deconfinement phase 
transition becomes explicit. 

It should be noted that the absence of real poles in 
the quark propagator is a sufficient but not a necessary 
criterion for confinement. In fact, in contrast to older 
studies |M]) the Schwinger-Dyson analysis of [M] seems 
to favor the existence of a real quark pole when a trunca¬ 
tion scheme beyond the rainbow-ladder approximation is 
used. It was also pointed out that the existence of com¬ 
plex conjugate poles in the nNJL model leads to thermo¬ 
dynamic instabilities [39] and only inhibits the decay of 
mesons into quark-antiquark pairs if additional prescrip¬ 
tions about the integration contour are made [53] . 

Nevertheless these models provide a nice and relatively 
simple test bed to investigate whether the pole structure 
of the quark propagator (“confining” or “non-confining” 
in the above sense) has a qualitative effect on the be¬ 
havior of Si. In this article, we therefore explore the 
behavior of Si in the NJL and nNJL models such that, 
with the appropriate choice of parameters, confinement 
is built-in or not in the models. Thereby we address the 
validity of Si as an order parameter for the confinement- 
deconfinement transition in these models. 

The remainder of this article is organized as follows: 
In Sec. El we introduce the NJL and nNJL models and 
explain how the chiral condensate and Si are obtained 
within these models. The corresponding results are pre¬ 
sented in Sec. [ml Finally, in Sec. m we draw our con¬ 
clusions. 


II. THE MODELS 


A. Local NJL model with infrared cutoff 

The Nambu-Jona-Lasinio (NJL) model [T71 [13] is a 
model of self-interacting fermions and was introduced in 


^ An exception is the propagator in m, which has no pole in the 
complex momentum plane but complex conjugate cuts. 
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the early 1960’s to describe the mass of nucleons through 
spontaneous chiral symmetry breaking. After the advent 
of QCD it was reinterpreted as an effective model for 
quarks, which acquire a constituent mass by the same 
mechanism (see Refs. for reviews). 

Throughout this article we work in Euclidean space, 
following the conventions of Ref. [U]. The NJL model is 
then given by the Lagrangian 

^NJL = + m)ip - Y f j, (3) 


where ip denotes a quark field with Nf = 2 flavor and 
Nc = 3 color degrees of freedom, and with bare mass m. 
The operator 0 is defined as 




(4) 


with 74 = 17 ° and the imaginary time variable X 4 = 
ix^ = it. The quarks interact by local four-point ver¬ 
tices, proportional to the coupling constant G. The 
interaction is invariant under chiral SU(2)l x SU{ 2 )fi 
transformations and consists of a scalar-isoscalar and a 
pseudoscalar-isovector part, where r denotes the triplet 
of the Pauli matrices in isospin space. 

In vacuum chiral symmetry, which is explicitly broken 
by the bare mass m, is also broken spontaneously through 
the self-interactions. In one-loop approximation (being 
equivalent to the mean-field or Hartree approximation), 
this gives rise to the dressed Euclidean quark propagator 


S{q) 


-^ + M 


(5) 


where M denotes the constituent quark mass, given by 
the gap equation 


M = m + G j ( 6 ) 

Here the trace is to be taken over the internal quark de¬ 
grees of freedom, i.e., Dirac, color and flavor components. 
One finds 


M = m + 8 NcG 


d^q M 
(27r)4 q^ + 


(7) 


Since the constituent quark mass also enters the right- 
hand side, the equation must be solved self-consistently, 
underlining its non-perturbative nature. 

The quark condensate is generally given by the expres¬ 
sion 

{ijjiP) = - J-^^ti{Siq)). ( 8 ) 

If we compare this with Eq. ([^ we find the simple relation 
/ T , 1 M — m 

{ipip) = - - —. (9) 


Applying Matsubara formalism, the above expressions 
can straightforwardly be generalized to nonzero temper¬ 
ature by the substitutions 



(74 —>■ Wra , (II) 


where a;„ = ( 2 n-|-l) 7 rr are fermionic Matsubara frequen¬ 
cies. Furthermore, recalling that the latter are a conse¬ 
quence of the antiperiodic boundary conditions fermions 
must obey in the imaginary time direction, the twisted 
boundary conditions defined in Eq. Q are easily imple¬ 
mented by the shift 


OJn 



( 12 ) 


So far we have ignored the fact that the integrals in 
Eqs. ([^ - (|^ and their extensions to nonzero tempera¬ 
ture and twisted boundary conditions diverge in the ul¬ 
traviolet and therefore have to be regularized. Since the 
NJL model is not renormalizable, the divergences can¬ 
not be absorbed in a redefinition of the parameters in 
the Lagrangian, so that the regularization scheme and 
the corresponding cut-off parameters must be viewed as 
a part of the model. 

As already mentioned, a second shortcoming of the 
NJL model in its original form is the fact that the model 
is not confining. Formally this can be seen most easily 
from the dressed propagator, Eq. ([^, which takes the 
form of the propagator of a non-interacting particle. In 
particular it has a pole in the time-like region at q^ = 
—M^, corresponding to the mass shell of a real particle. 

It is possible, however, to circumvent this problem by 
choosing a regularization scheme which avoids the poles 
of the propagator and, hence, the appearance of quark- 
production thresholds [301 El] • The basic idea is to re¬ 
strict the distance the confined quark can propagate by 
introducing an infrared cut-off, in addition to the one 
needed to regularize the UV divergences. The proper¬ 
time (PT) regularization turned out to be best suited for 
this task. To this end we use the identity 


1 


q^ + M^ 


dr e 


-T{q^+M^) 


(13) 


to write the propagator as an integral. If we now restrict 
the domain of the integration by introducing an IR cut¬ 
off and a UV one we find 




' WTW ■ 


, (14) 


where the original pole at is canceled by the 

numerator. Hence, we have removed the singularities 
from the propagator. As discussed before, this can be 
interpreted as confinement. 
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Applying this scheme to Eq. Q, the gap equation in 
vacuum becomes [30] 


M = m+^GM^ (r(-l, MVi) - r(-l, M\l)), 

(15) 

where r(a;,a;) is the incomplete gamma function. At 
nonzero temperature one gets 


2Nc 

M = m- ^GMT\^ 

7^3/2 Z—^ 






+y/TTMn[eiic{MnTuv) “ erfc(M„Tir) 


, (16) 


where = a/+ Cli^ and erfc(x) is the complemen¬ 
tary error function. From this we obtain the gap equa¬ 
tion for the mass with twisted boundary conditions if 
we replace by the shifted Matsubara frequencies 
given in Eq. (12). 


Having solved the gap equations, the chiral condensate 
is immediately obtained from Eq. Q, which also holds 
at nonzero temperature. In the same way we get the 
generalized condensate (' 0 ' 0)0 as 


= - 


G 


(17) 


from the solution of the gap equation with twisted 
boundary conditions. Si is then easily calculated from 
Eq. with n = 1. 

Finally, we would like to point out that the cancellation 
of the propagator poles in Eq. (14) is independent of the 


constituent quark mass. In particular, the quarks remain 
“confined” in the chirally restored phase where M = Oj^ 


B. Nonlocal NJL model 

The nonlocal Nambu-Jona-Lasinio (nNJL) model is de¬ 
scribed by the Euclidean Lagrangian 

G 

>C„NJL = + Tn)ip - —jaja, (18) 

where tp, m, and G are again the quark field, its bare 
mass, and a coupling constant, respectively. The nonlo¬ 
cality of the model is encoded within the nonlocal cur¬ 
rents ja{x) defined as 

ja{x) = J d'^yd'^z r{y - x)r{z - x)tp{y)Ta'ipiz), (19) 


^ In order to “cure” this problem the authors of Ref. m intro¬ 
duced a temperature dependent IR cut-off Tir(T), which diverges 
in the chirally restored phase, so that chiral and deconfinement 
phase transition coincide. 


with operators Tq G again corresponding to 

the scalar-isoscalar and pseudoscalar-isovector channels. 
The function r(x) is a regulator, which will be specified 
later. 


In the local limit, r[x) = 5{x)^ the integrals in Eq. (19) 
become trivial and we recover the standard NJL model, 
Eq. @. In general, however, the results get modified 
by the nonlocality. Most important, the one-loop quark 
self-energy, which for local interactions is constant in mo¬ 
mentum space, now becomes a momentum dependent 
function. As a result the dressed propagator takes the 
form 




( 20 ) 


where the function E(< 7 ^) replaces the constituent mass 
M in Eq. ([5|. 

For general nonlocal interactions the functional form 
of E( 9 ^) must be found by self-consistently solving a 
Schwinger-Dyson equation. In the present model, how¬ 
ever, a great simplihcation comes about from the fact 
that the interaction is separable, meaning that the four- 
point vertices in momentum space are essentially propor¬ 
tional to the product of the Fourier-transformed regula¬ 
tor functions. The form of the function 'S{q'^) is then 
simply given by 

=m + dr'^iq^), ( 21 ) 


where r{q^) is the regulator function in momentum space 
and (T is a constant, satisfying the gap equation 


d = G 


d'^p 


r 2 (p^)tr(S'(p)). 


Inserting Eqs. (20) and (21) this becomes 
d'^p 


d = 8NrG 


(27 




m + ar^{p^) 


+ (m + ar‘^{p^)y ’ 


( 22 ) 


(23) 


which, for a given r{q'^), must be solved for d. 

As before, the chiral condensate is given by Eq. (|^. 
However, if we compare this equation with the gap equa¬ 
tion (22) we see that the latter contains extra regulator 
functions in the integrand, so that in contrast to Eq. ([^ 
there is no simple relation between (ip'ip) and d (except 
for r(g^) = const., corresponding to a local interaction). 
Explicitly we have 


(' 0 '!/') = - 8 Vc 


d'^p m-\-(jr^{p^) 

(27r)‘* _|_ [jfi _|_ ^j.2(^p2y2 


(24) 


Finite temperature effects as well as twisted boundary 
conditions can again be incorporated through the substi¬ 
tutions Eqs. ( |I0| - (12) in the expressions above. 

If the function r(p) is chosen to drop off sufficiently 
fast at high values of p ^, the integral in Eq. (22) con¬ 


verges, so that no further regularization of the model is 
needed. Nevertheless, for the quark condensate, Eq. (24), 

























5 


there remains a quadratic divergence due to the bare 
quark mass. This problem also exists in QCD and re¬ 
quires a proper mass renormalization, see, e.g., Ref. m- 
In order to avoid such complications, we will restrict our 
numerical studies to the chiral limit, m = 0, where the 
quark condensate is convergent. 

Unlike the local NJL model, the dressed propagator 
in the nNJL model has a nontrivial pole structure. De¬ 
pending on the model parameters and in particular the 
regulator function, the pole positions given 

by 





= 0 


(25) 


can be real or complex. In the latter case, assuming 
that r{q^) is real for real the poles appear in complex 
conjugate pairs, which may be parametrized as 


= M^± iMT. (26) 

Following a quasiparticle picture, M can be interpreted 
as the constituent mass of the quark and F as its decay 
width [MlETlIle]. In this manner, a complex pole cor¬ 
responds to an unstable quasiparticle state, which could 
be interpreted as a manifestation of confinement]^ 

An important aspect in this context is that the pole 
structure may change as a function of temperature from 
“confining” to “not confining”. Hence, in contrast to the 
FT regularized NJL model, where the quarks are always 
confined, the nNJL model allows us to study not only the 
chiral phase transition, but also the deconfinement phase 
transition. 


consider the chiral limit, m = 0, and adopt the parame¬ 
ters of Ref. [47] 1^ 

G = 1.275 • 10"^ MeV-2 , (27) 

T„ = (240 MeV)-i , (28) 

Tuv = (905 MeV)-i , (29) 

which have been fitted to vacuum properties in the pion 
and rho-meson sector. With these parameters the vac¬ 
uum values of the constituent quark mass and of the chi¬ 
ral condensate per flavor are given by M = 358 MeV and 
{uu)^/^ = {dd)^/^ = —243 MeV, respectively. 



FIG. 1. Absolute value of the quark condensate ('i/’i/’) (solid 
line) and dnal qnark condensate Ei (dashed line) as functions 
of the temperature, calculated within the PT regularized NJL 
model in the chiral limit. The condensates have been normal¬ 
ized by the absolute value of {ijiip) at T = 0. 


III. RESULTS 

A. Local NJL model with infrared cutoff 


We are now ready to present our results for the behav¬ 
ior of the chiral condensate and the dual quark conden¬ 
sate as functions of the temperature. We begin with the 
PT regularized NJL model, introduced in Sec. |II A[ We 


^ This interpretation is probably too naive in various aspects. 
First, the underlying assumption is that the propagator can be 
Wick-rotated in the usual way, so that a pole at in 

Euclidean space corresponds to a pole at in Minkowski 

space. However, this property is spoiled by the existence of the 
complex conjugate poles themselves. Second, physical resonances 
have a cut along the real axis and complex poles are allowed 
only on the second Riemann sheet. One may, however, turn this 
argument around: Complex conjugate poles on the first Riemann 
sheet violate microcausality and positivity. Therefore such poles 
do not correspond to physical particle states, which in turn could 
be interpreted as confinement. 


At nonzero temperature, we obtain the results dis¬ 
played in Fig. The solid and the dashed line indicate 
the absolute value of {'ipip) and Ei, respectively, both 
normalized by the absolute value of ('i/'V’) at T = 0. 
We find that (t/’t/') continuously goes to zero at a crit¬ 
ical temperature Tc = 216.2 MeV, i.e., at this point a 
second-order chiral phase transition takes place. Ei, on 
the other hand, is very small at low temperatures and 
then smoothly rises, with a maximum slope at T = Tc- 
In QCD this would indicate a (rather broad) crossover 
from the confined to the deconfined phase. In the present 
model, however, the quarks remain confined, as discussed 
above. Hence, the rise of Ei has nothing to do with a 
deconfinement transition in this case, but is simply trig¬ 
gered by the change of the chiral properties. 


® The model of Ref. 1471 uses a local four-point interaction with 
the quantum numbers of a heavy-gluon exchange, parametrized 
by an effective coupling parameter tuq = 132 MeV. In the 
scalar-pseudoscalar channel this is Fierz equivalent to the La- 
grangian (j^ with G = 
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B. Nonlocal NJL model 


For our investigations within the nNJL model we take 
a Gaussian regulator 

(30) 


which drops off exponentially at large Euclidean mo¬ 
menta, so that no further regularization is necessary. In¬ 
serting this into Eqs. (21) and (25) one finds that the 
propagator has an infinite number of poles in the com¬ 
plex plane [44) . However, depending on the model 
parameters, there may be poles at real q^ as well. In the 
chiral limit, which will be considered in the following, the 
propagator has two real poles if the gap parameter a is 
smaller than a critical value (Ml 


^crit /7W ’ 

\l le 

while there are no real poles if ct > Ucrit ■ As discussed be¬ 
fore, the latter case can be interpreted as a manifestation 
of confinement. 

In our studies we therefore consider two qualitatively 
different sets of parameters, which are listed in Table 
m together with the corresponding value of Gcrit and the 
vacuum solution tto of the gap parameter. Erom these one 
can see that for parameter set A the quarks in vacuum 
are confined (in the above sense), while for set B they are 
not. As temperature increases, a decreases and finally 
vanishes at the chiral phase transition. Hence, for set 
A, there is also a deconfinement phase transition where 
a drops below Ocrit- For set B, on the other hand, the 
system is always in the deconfined phase. We are thus 
led to the question whether this qualitatively different 
behavior is also seen in the dual quark condensate. 


Set 

A [MeV] 

G [MeV^] 

UcHi [MeV] 

CTO [MeV] 

A 

760 

3.6- 10“® 

326 

404 

B 

914 

2.1 ■ 10“® 

392 

325 


TABLE I. Two sets of model parameters (regulator scale A 
and coupling constant G), the corresponding critical gap pa¬ 
rameter (tcrft, and the solution ao of the gap equation (22) 
at zero temperature. The calculations are performed in the 
chiral limit, m — 0. 


Our results for parameter set A are displayed in Fig.[^ 
Besides Ei and the chiral condensate we also show the 
mass M and the decay width F, according to the defini¬ 
tion in Eq. (26). Here we focus on the propagator pole 


with the lowest mass, since it contributes most signifi¬ 
cantly to the thermodynamics of the model. We can see 
that the pole, which is complex at low temperatures be¬ 
comes real at T =110 MeV, which, according to the in¬ 
terpretation discussed above, should be identified with 
the deconfinement temperature T^. The chiral phase 
transition, signaled by the vanishing of the mass and the 
chiral condensate, takes place in the same regime but at 



FIG. 2. Various quantities calculated with parameter set A as 
functions of the temperature: Mass (dotted line) and decay 
width (dash-dotted line), according to Eq. (261 applied to the 
first propagator pole; {4!’tp) (solid line), and Ei (dashed line). 


a slightly higher temperature, = 126.5 MeV. This is 
of course expected because d’(r) first drops below acrit 
aX T = Td before it vanishes completely at T = Tc- 
El is again very small at low temperatures and then 
rises significantly. Like in the PT regularized local NJL 
model, the maximum slope is found at the chiral transi¬ 
tion temperature Tc- However, since Ei rises smoothly, 
a relation to the deconfinement temperature Td cannot 
be totally excluded from the figure. 



FIG. 3. The same as Fig.[^ but for parameter set B. The de¬ 
cay width corresponding to the first propagator pole vanishes 
for all temperatures. 


We therefore turn to the results for parameter set B, 
which are shown in Fig. Just as for parameter set 
A, the rise of Ei occurs around the chiral phase transi¬ 
tion, signaled by the vanishing of the chiral condensate. 
However, in this case there is no deconfinement phase 
transition occurring. Hence, the order parameter like be- 
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havior of the dual quark condensate can only be related 
to the chiral symmetry and should not be regarded as a 
sign of confinement in the model. 

IV. CONCLUSIONS 


In the present article we have studied the dual quark 
condensate Si in local and non-local variants of the NJL 
model. In QCD, Eiis an order parameter for confine¬ 
ment [S], which is related to the breaking of the Z(3) 
center symmetry and as such equivalent to the standard 
“thin” PL in the limit of infinite quark masses. In the 
NJL model, on the other hand, there is no center sym¬ 
metry and no confinement, but nevertheless Ei behaves 
in a qualitatively similar way as in QCD (TB]. In Ref. [20] 
this was shown within a Ginzburg-Landau type analysis 
to be an effect of the chiral phase transition. 

Like the original NJL model, the models studied in 
the present paper do not have gauge helds and, hence, 
no center symmetry. However, there exists the possi¬ 
bility to have quark propagators without poles on the 
real Euclidean axis, which is often interpreted as a 
realization of confinement. The aim of our analysis was 
therefore to investigate whether the pole structure of the 
quark propagator, conhning or non-confining in the above 
sense, leaves imprints on Si. 

Specifically, we have considered three examples which 
all feature a chiral phase transition but have rather differ¬ 
ent confinement properties: a PT regularized local NJL 
model with infrared cutoff, which is always confining, a 
nonlocal NJL model with a deconfinement phase transi¬ 
tion at finite temperature, and a different parametriza- 
tion of the same model where the quarks are always de- 
confined. 

We find that Si behaves almost identically in all three 
cases, namely it rises most steeply at the chiral phase 
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transition temperature, just like in the NJL model. In 
particular, we do not see any effect related to a change 
of the conhning properties of the propagator. Although 
in one parametrization of the nNJL model there is a de- 
conhnement transition and the rise of Si roughly falls in 
that region, this must be seen as a coincidence because 
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if this is also true in QCD, it could explain the approxi¬ 
mate coincidence of chiral and deconhnement crossovers, 
observed in lattice calculations. 
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